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Modern intense ultrafast pulsed lasers generate an electric field of sufficient 
strength to permit tunnel ionization of the valence electrons in atoms^ This process 
is usually treated as a rapid succession of isolated events, in which the states of the 
remaining electrons are neglected^. Such electronic interactions are predicted to be 
weak, the exception being recoUision excitation and ionization caused by linearly- 



polarized radiation^. In contrast, it has recently been suggested that intense field 
ionization may be accompanied by a two-stage 'shake-up' reaction"*. Here we 
report a unique combination of experimental techniques^ " * that enables us to 
accurately measure the tunnel ionization probability for argon exposed to 50 
femtosecond laser pulses. Most significantly for the current study, this 
measurement is independent of the optical focal geometry''*, equivalent to a 
homogenous electric field. Furthermore, circularly-polarized radiation negates 
recoUision. The present measurements indicate that tunnel ionization results in 
simultaneous excitation of one or more remaining electrons through shake-up'. 
From an atomic physics standpoint, it may be possible to induce ionization from 
specific states, and will influence the development of coherent attosecond XUV 
radiation sources^". Such pulses have vital scientific and economic potential in 
areas such as high-resolution imaging of in-vivo cells and nanoscale XUV 
lithography. 

The Noble gas atoms are a natural choice for studying electron dynamics in an 
ultrafast intense laser pulse, and a number of pivotal experimental studies have 
established this rich and vibrant field of research, revealing much about the complex 
nature of photoionization^'. However, the majority of these studies have employed 
linearly polarized radiation, and used the predictions of tunnelling theory^ (generally the 

12 

ADK treatment ) in which the target is confined to the ground state. The current 
investigation was stimulated by the theoretical work of Eichmann et al^^ and Zon et af, 
and reveals experimental evidence for a hereto unobserved accompanying excitation 
mechanism. 

hi the ultrafast regime, the optical pulse duration is of the order of femtoseconds 
(1 femtosecond = 10"^^ s) and ionization proceeds either by a multiphoton perturbative 
process or a nonperturbative strong -field process, generally accepted as being 



described by tunnel theory: for a recent review see Popov . The present work is 
concerned purely with strong field processes (peak intensity greater than approximately 
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10 TW/cm , where 1 TW =10 W) where the rate of tunnel ionization is governed by 
the frequency and strength of the radiation-induced electric field and the binding energy 
and quantum state of the ion and electron(s). 

Immediately following ionization, the 'drive pulse' electron is in a Volkov state 
and initially is fully correlated with the parent ion. As the electron is accelerated by the 
laser field, the electron trajectory is determined by the ellipticity of the field and the 
phase at which ionization occurred. Electron impact excitation or subsequent 
ionization can arise in a linearly polarized laser field through the intriguing 
phenomenon of recoUision. RecoUision is the key mechanism for attosecond XUV pulse 
generation, as the kinetic energy of the electron is dissipated photonically if the electron 
is recaptured by the parent ion^^. However, in the present work, we make recoUision 
events negligible by employing circularly polarized light: the absorption of a large 
number of photons transfers considerable angular momentum to the liberated electron, 
preventing it from returning to the ionic core ' , thus the masking effect of recoUisional 
excitation and ionization are negated. 

While strong-field ionization and recoUision in ultrafast laser pulses is well 
documented, minimal theoretical and essentially no experimental studies have 
investigated the possibility of simultaneous excitation of the parent ion during tunnel 
ionization. In the case of multiphoton ionization, resonant excitation processes have 
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been identified " , however in the present work, the intensity of the laser field is 
sufficient that tunnel ionization occurs with a dramatically enhanced probability. The 
contemporary work of Zon"^'^"*'^^ invoked the idea of 'inelastic tunnelling' whereby the 
parent ion is left in an excited state following the ionization of one of N identical 
valence electrons. Excitation is through 'shake-up', first employed by Carlson^ to 



explain single UV photon absorption leading to the ionization of a first and the 
excitation of a second electron. The ionization event diabatically distorts the electron 
wavefunctions, resulting in the excitation of a bound second electron. Using this 
concept, Zon et af'^^ and Komev et af^ have derived a general expression for the rate 
of tunnel ionization (TI) of an atom with simultaneous excitation of the lowest lying 
ionic states. This treatment relies on the sudden approximation (SA), which is valid 
provided accessing the excited states in Ar*^'''^'^'^ requires considerably less energy than 
the ionization potential of Ar'''^. Komev et af'^ show that while this has negligible 
influence on the ionization probability for the Ar'''^ — > Ar*^''^^^^ process, the creation of 
excited states in the product Ar^'^^'^^ ion strongly influence the Ar^'^^'^^ Ar*-''^^^^ 
ionization probability. The consequence of this mechanism is illustrated in figure 1. 
Furthermore, by allowing for all combinations of excitation and ionization, the 
significant role of simultaneous excitation and ionization during tunnelling in the laser 
field is highUghted. To distinguish from standard sequential tunnel ionization (TI), we 
refer to such processes as multi-electron tunnel ionization (METI). 

The analogous process of 'shake-off resulting in further ionization rather than 
discrete excitation was originally proposed a decade ago by Fittinghoff et al in an 
attempt to quantify the many-orders of magnitude discrepancy between the predicted 
sequential TI jdeld and the observed He ion yield by laser radiation . It has been 
through careful experiments emplo5dng linear and circular polarization (see for example 
the workofGuoe^a/ ) and the COLTRIMS technique ' that recoUision has finally 
been accepted as the mechanism responsible. More recently, the shake-off ionization of 
helium exposed to ultrafast laser pulses has been theoretically investigated using the S- 
matrix technique. Becker and Faisal calculated the ratio of the rates of shake-off to 
recoUision, and demonstrated that shake-off was many orders of magnitude weaker^^. 
As a result, the presence of shake-off in the strong-field ionization of atoms has been 
generally dismissed for the case of linear polarization. In the present work, we propose 



that the relevance of shake-up excitation be re-examined in the case of argon, given the 
validity of the SA is significantly different from the case of heUum discussed above. 
The far lower-lying excited states in Ar'^^ (q = 1 to 6) are far more readily accessed 
during TI, in contrast to the extremely high excitation energies of He^, irreconcilable 
with the SA. 

In general, previous experimental measurements of ion yield as a function of laser 
intensity are a convolution of the ion signal with the focal volume producing the signal. 
By simply changing the energy of the laser pulse, the spatial distribution of laser 
intensity also changes. Frequently, the complexity of this situation is compounded by 
diffraction associated with the spatial profile of the laser ; a direct comparison with 
theory is thus impossible without introducing the specific experimental geometry. A 
novel solution to this 'volume variation' problem, referred to as intensity selective 
scanning (ISS), was demonsfrated by Van Woerkum^; the relevance of ISS is described 
in Methods, and is illustrated in figure 2. The ion yield recorded at each position Zf, the 
distance of the spectrometer axis from the focus, still inherently depends on the intensity 
distribution in the {x, y) plane. Two numerical methods have been proposed for 
removing spatial integration over each 'slice' through the focal volume: Walker et af" 
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described a deconvolution routine, and Goodworth et al proposed a tomographic 
method. Both recover focal geomefry-independent probabilities of ionization; the 
authors have extended the former to experimentally realistic non-Gaussian focal 
conditions . 

In figure 2, we present a measurement of the spatial distribution of ions generated 
when argon is exposed to circularly polarized 50 fs pulses over the intensity range 100 
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TW/cm to 100 PW/cm . This experiment was performed using the ASTRA Laser 
Facility at the Rutherford Appleton Laboratory, as detailed in Methods. The Ar^^ = 1 
to 6) ion yield was measured as the laser focus was translated by a computer-driven 



motion stage, exposing tlie detector to different 'slices' through the focal volume: the 
spatial selectivity of the time-of-flight mass spectrometer (TOFMS) (grey arrow) is 
indicated by the narrow slice in the (x,y) plane. 

The partial probability of ionization PPI to Ar*^ = 1 to 6) shown in figure 3 is 
recovered from the ISS data (figure 2) using the technique pioneered by the authors^. 
Importantly, we totally remove any instrument dependence, including the enhancement 
of detector gain with charge state. The distinction 'partial' is necessary as the 
deconvolution method is only valid at intensities below saturation (as indicated by the 
dashed line), irrespective of diffi-action effects. The saturation intensity is defined as 
described in Methods. The PPI for each ionization state (figure 3) can be directly recast 
in terms of the universal conserved probability of ionization (CPI), which embraces all 
charge states. The CPI for each charge state is presented in figure 4, and is discussed in 
Methods. 

The spatially homogenous intensity-dependent predictions of Kornev et al are 
perfect for direct comparison with the measured CPI data. Such a comparison is 
presented in figure 4(a)-(f), where the data points are the present measurements, the thin 
lines the calculated probabilities of producing Ar''"^ from the Ar*''"'^'^ ground state only 
(sequential TI), and the thick lines are the METI predictions including shake-up. Note 
the experimental CPI has been normalized to the theoretical predictions via the rising 
edge of the Ar^ CPI; this first ionization stage will not be influenced by excitation. 
Figure 4 clearly highlights the significant difference between theory and experiment 
when only sequential TI from the ground state is considered (thin line), lying generally 
well outside the experimental uncertainty. This is illustrated in figure 4 by 
representative error bars. However, when shake-up is included in the METI prediction, 
an excellent agreement is found for all charge states as indicated by the thick lines in 
figure 4. While other possible excitation mechanisms could also be considered, for 



example the 'way-out' excitation, whereby the departing electron scatters from another 
bound electron , the observed agreement is strongly supportive of the shake-up model. 
This benefits from the repeated diabatic distortion of the valence electronic 
wavefunctions during tunnel ionization at high intensities over a number of optical 
cycles at intensities greater than 1 PW/cm . Thus is an important contrast to single 
photon shake-off , as the cumulative influence of the oscillatory laser field need be 
considered. 

To conclude, we have recorded for the first time strong evidence for the presence 
of considerable atomic excitation during tunnel ionization by a 790 nm 50 fs circularly 
polarized laser pulse focused to intensities in excess of lOOTW/cm^ Such an 
observation is made without the need to embrace recoUision processes. The significance 
of the agreement between our experimental observations and recent theoretical 
predictions indicate that excitation during ionization must be considered irrespective of 
recoUision processes, referred to as multi-electron tunnel ionization. Such excitation 
will also occur in a linearly polarized laser field, and is expected to be even more 
influential in a 5 fs few-cycle laser pulse . As the proposed method for generating 
intense bursts of XUV radiation relies on accurately controlling high harmonic 
generation in few-cycle intense laser fields excitation during tunnelling has a major 
bearing on the emerging field of optical attosecond physics. Accurate manipulation of 
the liberated electron motion during the infrared 'drive pulse' and the occupation of 
initial and transient electronic states will allow tuning of the spectrum of the XUV 
pulse. 

Methods 

Intensity Selective Scanning. We are interested in ionization and excitation processes 
that saturate at intensities less than 100 PW/cm , therefore the 30 mJ 790 nm 50 fs laser 



pulses generated by the ASTRA Laser Facility (UK) need only be softly focussed (f/ll 
optics) to generate a peak intensity in excess of 100 PW/cm". Indeed, multiple 
ionization occurs over tens of millimetres. By measuring the Ar*^ = 1 to 6) ion yield 
with a tightly apertured (250 microns) ion time-of-flight mass spectrometer, only those 
ions generated within a narrow spatial (and therefore intensity) window are detected. 
Then, as the focusing optic is translated, the spectrometer is only exposed to ions 
generated by a well-defined laser intensity " ' , as illustrated in figure 2. Throughout 
this measurement, the argon gas pressure is low enough so as to avoid space-charge 
effects, tested by repeating the ISS measurements as a function of target gas pressure. 
The narrowing of the spatial distribution of ions as a function of Zf for increasing n is 
representative of the increasing ionization potential, where the highest ionization state 
presented (Ar^^) is peaked at Zf = mm where the intensity is at a maximum. 

Deconvolution of a non-Gaussian focus. The ion yield measured by intensity selective 
scanning depends very specifically on the spatial distribution of intensity within our 
focal volume. To make this measurement more universal, we remove this dependence 
through a deconvolution technique , requiring the measured ion yield and theoretical 
on-axis intensity as a function of focal position as inputs. The resulting partial 
probability of ionization (PPI) is independent of the variation of the signal -producing 
volume. We are also able to account for the unavoidable diffraction of the laser pulse 
by solving the Huygens-Fresnel diffraction integral. The PPI is equivalent to the 
response of a single atom to a spatially infinite laser focus, directly comparable to 
theoretical predictions, which tend to be presented in terms of a homogenous laser field. 
The PPI results are valid up to saturation, however at higher intensities the 
deconvolution breaks down, and the PPI is defined as unity. By measuring the first 
derivative of the PPI, the saturation intensity is defined as the point at which the 
gradient falls below the equivalent uncertainty on the PPI. The uncertainty is estimated 
by calculating the statistical deviation in the recorded average time-of-fiight signal, and 



the upper and lower confidence interval corresponding to 1 s.d. was propagated through 
the deconvolution procedure. Furthermore, the probability at which saturation occurs 
depends on the quantum efficiency of the detector, which varies with charge state. The 
deconvolution method can be used to remove this dependence by normalising the 
saturation PPI to unity. 

Conservation of probability. Consider the PPI(n) for « = 1, 2 as shown in figure 3. At 
low intensity, PPI(l) is small, and increases with intensity up to saturation at an 
intensity of 1 PW/cm . However, as is often the case in atomic ionization, PPI(2) is 
nonzero at this intensity. The condition for conserving probability is that the sum of 
probabilities is less than unity below the saturation intensity of PPI(l) or equal to unity 
above the saturation intensity of PPI(l), thus at an intensity greater than ~1 PW/cm , the 
conserved probability of ionization to Ar^, CPI(l) must be less than unity as PPI(2) is 
nonzero. This definition is extended to an N electron system in the present work. The 
uncertainty in the CPI is estimated by propagating the uncertainty in the PPI through the 
deconvolution procedure. 
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Figure 1 Tunnel ionization and resulting excitation mechanisms in the Coulomb 
potential of argon. Regular sequential tunnel ionization: a. at an intensity of 1 
PW/cm^ the ground state CP312) of Ar"" has a significant probability to tunnel 
ionize (1). b. As the intensity is increased to 3 PW/cm^, the ground state (^P2) of 
Ar^"^ can ionize (2). Considering multi-electron processes, the tunnel ionization 
is significantly influenced by excitation, c. As the ^P3/2 electron tunnels from the 
Ar"" (1), a significant population is transferred from the ^P2 to ^Pi excited state in 
the Ar^"" ion (2). d. A lower intensity (2 PW/cm^) is then required to ionize 
population in the ^Pi state (3), while the ^P2 state ionizes at 3 PW/cm^ (b). The 
generation of Ar^"" can then further excite the Ar^"" ion (4). We measure the 
dependence of the probability of ionization on the laser intensity and charge 
state. 




Figure 2 Illustration of the intensity selective scanning (ISS) technique^ "® and 
measured ion yield data. As the laser focus (bottom) is translated with respect 
to the narrow detection aperture of the TOFMS (slot perpendicular to Zf, the 
distance from the focus), the ion yield for all charge states of Ar'^'" {q = 1 to 6) is 
measured. The laser pulse propagates in the z direction, and the focus is 
translated parallel to this axis. 
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Figure 3 Partial probability of ionization (PPI) to Ar'^'^ {q = 1 to 6) as a function of 
spatially independent laser intensity. As circularly polarized radiation is 
employed, this data is free from the influence of recollision effects. The partial 
probability is recovered from the ISS data in figure 2 using the method 
described in reference [7]. The uncertainty is estimated by calculating the 
statistical deviation in the recorded average time-of-f light signal, and the upper 
and lower confidence interval corresponding to 1 s.d. was propagated through 
the deconvolution procedure. The dashed line indicates a PPi of unity, whereby 
ionization is saturated. 
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Figure 4 Conserved probability of ionization (CPI) to Ar'^'" {q = 1 to 6) as a 
function of spatially-independent laser intensity, (a) Ar"", (b) Ar^"", (c) Ar^"", (d) 
Ar"*"^, (e) Ar^"^ and (f) Ar^"". Data points: current results derived^ from the data in 
figure 3. Thin lines: theoretical predictions^^ of sequential multiple tunnel 
ionization (Tl) from ionic ground states only, equivalent to ADK theory^^, thick 
lines: full multi-electron tunnel ionization (MET!) treatment allowing for shake-up 



excitation during tunnelling. The uncertainty in the CP! is estimated by 
propagating the uncertainty in the PPI through the deconvolution procedure. 



